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Abstract 

We consider the simplest extension of the standard electroweak model by one sterile 
neutrino that allows for neutrino masses and mixing. We find that its leptonic sector 
contains much less free physical parameters than previously realized. In addition to the 
two neutrino masses, the lepton mixing matrix in charged current interactions involves 
(n—1) free physical mixing angles for n generations. The mixing matrix in neutral current 
interactions of neutrinos is completely fixed by the two masses. Both interactions conserve 
CP. We illustrate the phenomenological implications of the model by vacuum neutrino 
oscillations, tritium (3 decay and neutrinoless double (3 decay. It turns out that, due to 
the revealed specific structure in its mixing matrix, the model with any n generations 
cannot accommodate simultaneously the data by KamLAND, K2K and CHOOZ. 

PACS: 14.60.Pq, 14.60. St, 23.40.-s 

Keywords: lepton mixing, neutrino oscillation, beta decay 



liaoy@itp.uni-leipzig.de 



1 



1 Introduction 



The experiments with solar, atmospheric, reactor and accelerator neutrinos have now 
provided evidence that neutrinos have mass and mix [1, 2]. These experimental results 
can be best understood in terms of neutrino oscillations [3]. This implies unambiguously 
that the standard model (SM) of electroweak interactions has to be extended at least 
in the leptonic sector. Most phenomenological analyses performed so far assume three 
active massive neutrinos and a three by three unitary mixing matrix for charged current 
(CC) interactions of leptons. The leptonic sector then contains two differences of neutrino 
masses squared, three mixing angles plus three CP-violating phases. It has been shown to 
be capable of accommodating all neutrino data except the LSND result [4] , which however 
was not confirmed by other short baseline experiments [5] and remains to be clarified by 
MiniBooNE [6] in the future. From the model building point of view, such a scenario 
requires new degrees of freedom to be added to SM, either some heavy neutrinos [7], or 
a Higgs triplet which also develops a vacuum expectation value [8]. In the first case, the 
three by three mixing matrix amongst three light, active neutrinos is only approximately 
unitary to the extent that their mixing with heavy neutrinos can be ignored in the analysis 
of current data. Thus it amounts to a low energy effective theory of a fundamental one 
which could be much more complicated. In the second one, the existence of a non- 
doublet Higgs boson is always severely constrained by precision electroweak data and the 
null result of direct Higgs searches. 

There are also attempts to incorporating the LSND result by including explicitly a 
sterile neutrino into the mixing scheme which was introduced earlier in the other context 
[9]. But they are found to be disfavored by the experimental data either because of 
the tension between the positive result of LSND and the negative ones by other short 
baseline experiments or because of the rejection of sterile neutrino's involvement in solar 
and atmospheric data [10]. 

In this work, we put aside the LSND result as in most studies and ask whether it is 
possible to understand the neutrino data in a minimal extension of SM. The extension 
is minimal in the sense that it introduces the least numbers of new degrees of freedom 
and free physical parameters into a fundamental theory without endangering precision 
electroweak data. With this in mind, one possibility would be to introduce one sterile 
neutrino to the SM of three generations. Such kind of models were systematically studied 
long time ago in the pioneering work of Ref. [11]. They were also considered without in- 
cluding Majorana mass terms in Ref. [12]. According to the analysis in Ref. [11], this would 
introduce 5 mixing angles and 3 CP-violating phases, in addition to two neutrino masses 
(with the other two being massless). There should thus be some room to accommodate 
the mentioned neutrino data that essentially call for two mass-squared differences and 3 
independent mixing angles. However, as we shall analyse in the next section, there are 
actually only 2 mixing angles and no CP-violating phase in the lepton mixing matrix after 
we take into account completely the constraints from the texture zero that appears in the 
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original neutrino mass matrix [13]. Furthermore, as we shall show in section 3, due to 
the special structure in the mixing matrix exposed in section 2, it is even not possible to 
accommodate simultaneously the data by KamLAND, K2K and CHOOZ [14, 15, 16] for 
any number of n generations, although in that case there are (n — 1) free mixing angles 
at our disposal. The same structure also results in the vanishing leading contribution to 
neutrinoless double j3 decay for any n which in principle is allowed to occur due to the 
Majorana nature of the neutrinos. Our results are summarized and conclusions are made 
in the last section. 



2 Parametrization of mixing matrix 

In this section we first describe the leptonic sector of the n-generation SM extended by 
no sterile neutrinos. Then, we specialize to the simplest case of no = 1 to parametrize the 
lepton mixing matrix and count its independent physical parameters. 



2.1 The model 

The only new fields compared to SM are the Uq sterile neutrinos that we choose to be 
right-handed without loss of generality, s Rx , x = 1,---,uq. It is sufficient for us to 
concentrate on the leptonic sector of the model that contains the standard n generations 
of the doublets, F^ a = (riLa, fLa) T , and of the charged lepton singlets, f Ra , a = 1, • • • , n. 
Here L, R refer to the left- and right-handed projections of the fields in terms of P LjR = 
(1 =F 7 5 )/2. The kinetic and gauge interaction terms of the leptons are 

£k = nE^ipn La + SR^iPs Rx + fjftfa 



r 9 7 



^ 1 (!) 
nTaYn La + f a Y (-P L + 2s 2 w ) f a 



where summation over indices a, x is implied, g, e are respectively the couplings of SU (2) l 
and C/(1)em, an d cw = mw/mz with m wz being the masses of W ± , Z. 

Since the sterile neutrinos are neutral under SU(2) L x U(1) Y by definition, they can 
have bare mass terms of Majorana type, 

-C SR = ^M xy ~4~ x s Ry + \M: y s^s c Rx (2) 

where ip c = C"f°ip* stands for the charge-conjugate field of ip with C = i'j ^ 2 satisfying 
C = -C f = -C T = -C- 1 and CY T C = Y- We denote sg = (s R ) c for brevity. The 
no x no complex matrix M is symmetric due to anticommutativity of fermion fields, but 
is otherwise general. Together with the Yukawa terms of the leptons, 



= y f ab F La <pfRb + Vax F La(pSRx + h.C. (3) 
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where tp is the Higgs doublet field that develops a vacuum expectation value, (tp) = 
(0, l) T v/y/2, and tp = ia 2 {p*, the lepton mass terms become 



f L m f f R + nZDs R + h.c] + - [s c R Ms R + h.c] (4) 



where = and D = y n v/ \/2 are nxn and nxn complex matrices respectively. 

The charged lepton mass terms are diagonalized as usual by biunitary transformations, 

f L = X L £ L , f R = X R £ R , X^ R = X ] LR , (5) 

with 

X\m f X R = m e = diag(m e , m M , m T , ■ ■ •) (6) 

being real and positive. We shall denote the mass eigenstate fields of the charged leptons, 
£l,r, by the Greek indices, a, (3 — 1, 2, • • • , n, corresponding to the electron, muon, etc. 

The kinetic and mass terms for the neutral leptons contain the fields ul, s r and s^. 
To diagonalize, we first rewrite them uniformly in terms of the fields ul, s r and their 
charge-conjugates n£, s R . Using ip c i^^d^x C = —^^(xi^'tp) + X^7 M ^V^> wnere t ne total 
derivative term can be ignored in the Lagrangian, and ip c x C — X>Pi the terms become 



,c 



S R 



(7) 



where the (n + n ) dimensional, symmetric mass matrix in the new basis is 

with n being the zero matrix of n dimensions. For n > tiq, which covers our interested 
case n — 3, no = 1 later on, it contains a zero eigenvalue of degeneracy (n — no) and 2tiq 
eigenvalues which are non-zero and nondegenerate for general parameters D, M. Without 
changing the diagonal form of the kinetic terms, we make a unitary transformation 



n 



c 



L ) = Yu R , Y- 1 = Y\ (9) 
sr J 



which also fixes the transformation of the conjugate fields, 



n ,c = Y*" C r, (10) 

S R J 



such that 

Y T m n Y = m u = diag(0, • • • , 0, m n _„ 0+ i, • • • , m n+no ) (11) 
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with the nonvanishing masses being real and positive. We shall denote the mass eigenstate 
fields of the neutral leptons, u R , by the Latin indices, j,k = 1, 2, • • • , n + n . Then, 

a + C = \ (mfrR + v&hf) - \ {f R m v v R + mrrvg) (12) 
which may be put in the compact form 

CL + £L = \H$-rn v )v (13) 
by introducing the Majorana neutrino fields 

v = vr + Vr (14) 

satisfying v c = v. 

Now we express the interactions of leptons in terms of the fields with a definite mass. 
There are no changes in + ^em f° r the charged leptons. The CC interaction becomes 
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VgW^tz.rfvj + VftW^-flLft] (15) 
where the n x (n + n ) matrix 

n 

V S = E(*i)/*% (is) 



a=l 



is the leptonic analog of the mixing matrix Vq KM in the hadronic sector. It is important 
to notice that only the first n rows in Y are involved in the CC mixing matrix V c [11] 
because the remaining n rows are associated with the sterile neutrinos sr which do not 
enter any interactions. Due to this and unitarity of Xl and Y, we have 

V C V C ^ = l n , (17) 

but V°W C ^ 

ln+no- Actually, the latter appears in the neutral current (NC) interaction 
for the neutrinos 



-nc - x—VfiZ&'fPz.Vi (18) 
by the relation [11] 



V% = E YuStaYZ = E J2Y bk (X L ) ba (Xl) aa Y a ) = (V C V% (19) 

a, 6=1 a, 6=1 « 

Using tP c ^PlX C — —XI^Pr^^ v f — v ji an d Hermiticity of U N , the interaction can also 
be cast in the form, 

£ nc = -£-Z,TT kl " [i Im V$ - 7 5 Re V*\ »j (20) 
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In addition to the condition (17), V c satisfies a relation that will be important in 
its parametrization. The original neutrino mass matrix m n has an n dimensional zero 
submatrix in its left-upper corner which is protected by gauge symmetry of the model. 
Then, eq. (11) implies that 

(Y*m v Y^) ab = 0, a,b=l,---,n (21) 

Multiplying it by (X| / ) aa (X^) /3b , summing over a, b and using eq. (16) leads to the matrix 
relation [13] 

V c m u V CT = n (22) 

Note that the above holds irrespective of n > n when there are massless modes or n < n 
when there is none for general parameters D, M. 

2.2 The mixing matrix 

From now on, we shall restrict ourselves to the case of no = 1. We offer two ways to 
construct the mixing matrix V c , one by using the constraints (17,22) and the other 
by explicit diagonalization. This is then followed by counting the independent physical 
parameters contained in it. 

We start with the case of n = 1 in which the two neutrinos are generally massive and 
nondegenerate with masses denoted by < m_ < m + . The mixing matrix is a one-row 
complex matrix, V c = (ne ipl , r 2 e ip2 ). Eqs. (17,22) read 

r\ + r\ = 1, m_r\ + m + rle 2i{p2 - pl) = (23) 

whose solution gives 



and then, 



V — e pi (c m , iis m ), c m — a , s m — a (24) 

y m + + m_ y m + + m_ 

The global phase in V c can be absorbed into the electron field while the two signs can be 
interchanged by flipping the sign of the field v 2 - Both mixing matrices are thus uniquely 
fixed by the two neutrino masses, in contrast to the claim in Ref.[ll] that there are one 
free mixing angle and one free CP-violating phase. Note that the nonreality of V C ' N is not 
a sign of CP violation. If we want, we can absorb the i into the field z/ 2 without changing 
its mass term so that V C,N are real. But this will introduce a nontrivial "creation phase" 
into it which was unity before rephasing [17]. 

The above result may also be obtained by direct diagonalization of the mass matrix 
in eq. (8) where D and M are two complex numbers. The unitary matrix Y in eq. (11) 
may be parametrized generally as 
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with c 2 + s 2 = 1 and all parameters being real. Denoting D = \D\e tSl , M = \M\e tS2 , eq. 
(11) then yields 



m± 



e i2 7l = _ e . 



(27) 



'|M| 2 + 4|D| 2 ± |M| 

3«272 _ _[_^ e — -i(<5i — <5 2 ) e *27o _ _|_^ e -i<5i 

and c, s being identified with c m , s m in eq. (24). Since only the first row of Y appears 
in V c , the relevant phase is e^ 72-7 ^ = ±i independently of 5^2, and the preceding result 
is reproduced. 

For the case of n > n = 1, we generally have 2 massive modes with masses m_ < m + 
that we arrange to be the last two in the order of increasing mass and (n — 1) massless 



modes. Assuming (V ) 



T ■ P 
a >3 



'l>a 



' , eq. (22) implies that 



m.r a , n r p , n e i ^ n+p ^ + m + r a , n+1 r p , n+1 e i ^+ 1+p ^ = 



(28) 



Note that the constraints from the off-diagonal elements (a ^ f3) are not independent but 
just resolve the separate two-fold ambiguities in the constraints for the diagonal elements 
(a — (5) to an overall two-fold ambiguity: 



fj 7 ^, e ipa ' n+1 = ±ie ipa '' 



(29) 



with r m = . V c can then be factorized as 



= VU 



where V is an n x n matrix whose a-th row is 

(r -.P*' 5 ". 1 ... r , pVa.n-l r . /I I r „ipa,n\ 

y a,l i i 1 a,n— 1° ; ' a,n V ~ ' m / 

and C/ is an n x (n + 1) matrix 



U 



ln-l 0( n _i) x i 0(„_1) X 1 

Olx(n-l) C m i^Sm 



(30) 



(31) 



(32) 



with c m , s m again given in eq. (24). l/W = l n together with eq. (17) gives, VV^ = l n . 
Since V is a square matrix whose entries are otherwise arbitrary, it is a general unitary 
matrix of dimension n. The above structure yields 



V 



N 



L n-1 



~ J r%C m S m S 



(33) 



Thus the NC interactions of neutrinos do not contain any free mixing angles or CP- 
violating phases that may appear in their CC interactions. They conserve CP, and the 
off-diagonal interactions occur only between the two massive neutrinos. 
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This result can be confirmed by direct diagonalization of the mass matrix in eq. (8). In 
the first step, we decouple the (n — 1) massless modes. Now D is an n-row column matrix 
whose entries are parametrized as D a = \D a \e %5a . The magnitudes define a vector in lZ n , 
which is rotated to the n-th axis by the rotation R . Denoting E = diag(e~ l<51 , • • • , e~ tSn ) 
and the (n + 1) dimensional unitary matrix, 

Y =^ RoEo)T J (34) 

we have 

Y T m n Y =| \D\ j (35) 

\D\ M 

with \D\ = \ D a \ 2 - The problem has thus been reduced to the case of n = uq = 1 
that we treated earlier. The required two-dimensional unitary matrix, now denoted as yi, 
for diagonalizing the submatrix is given in eqs. (26, 27). Then, 

Y T m n Y = diag(0,---,0,m_,m+) (36) 

where Y — Y Y 1 and 

y, = ( m ) (37) 

Noting the block diagonal form of F ,i, eq. (16) gives the same factorized form of V° as 
in eq. (30) with the same U as in eq. (32) but now 

V = X|(i? E ) t e i(70+7l) (38) 

Since Xl is a general unitary matrix, so is V. 

Now we count the free physical parameters contained in V c . An n dimensional unitary 
matrix V may be parametrized as a product in any arbitrarily specified order of the n 
phase matrices, e a (u a ) (a — 1, • • • , n), and the n{n — l)/2 complex rotation matrices in 
the (cx,f3) plane, oJapiQaPiVap) (n > (3 > a > 1) [11]. Here e a (z) is obtained by replacing 
the a-th entry in l n by the phase z, and 

uiap(0a(3, = e a (e l ^)R a p(6 a p)e a (e- % ^) (39) 

where R a f3(9 a p) is the usual real rotation matrix through angle 9 a p in the (a, (3) plane. 
We choose the order of products in such a way that it fits our purpose here: 

V = e n {u n ) \l Van) (^ X 1 ^ =V (^ X 1 ^ , (40) 

where X is the general (n — 1) dimensional unitary matrix containing (n— l) 2 real param- 
eters and the u an factors are ordered from left to right in increasing a. From the revealed 
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structure in eq. (30), the matrix containing X can be pushed through the matrix U to 
be absorbed into the massless neutrino fields without causing any other changes in the 
Lagrangian. This leaves us with (2n — 1) real parameters in Vq. 

However, not all parameters in V are physical. For brevity, we denote ui an (9 an , ip an ) = 
u a , R an (6 a n) = Ra, e a (e ilpan ) = e a , and e Q (e~^ a ") = e*. For n = l,V = ei(iti) can be 
absorbed into the single charged lepton field leaving no free physical parameters in V° as 
we found earlier. For n = 2, we have 

V = e 2 (u 2 )e 1 R 1 (9 1 )e\ (41) 

where e 2 (u 2 )e\ can be absorbed into the two charged lepton fields while e\ can be pushed 
through U and absorbed into the massless neutrino field, leaving behind one physical 
mixing angle in this case. For n > 3, we combine the two adjacent u;'s as follows, 

LO a uj a+ i = e a R a e a e a+ iR a+ ie a+1 = e a e a+ iR a R a+ ie a e a+1 (42) 

because [e a ,ep\ = 0, and [e a ,-Rg] = for a ^ {3 and a < n. The sequence can be 
continued such that 





n-l 


n-l 


"n-l 




Gn(^n) C« 
a=l 








13=1 





Again, the right factor commutes with U to get absorbed by the (n — 1) massless neutrino 
fields and the left one by the n charged lepton fields. 

To summarize, the lepton mixing matrix V c contains (n — 1) free physical mixing 
angles and no CP-violating phases, and may be parametrized as in eq. (30) with U given 
in eq. (32) and 

n-l 

V = II Ran(0 a ) (44) 

Q=l 

The angles 9 a describe the mixing of the (n — 1) massless neutrinos with the subsystem 
of the two massive ones. The mixing matrix in NC interactions of neutrinos does not 
contain these free parameters but is fixed by the two neutrino masses. It is off-diagonal 
only in the massive subsystem and also conserves CP. It is interesting that the above 
counting for V c happens to be the same as the one in Ref. [12] where the bare Majorana 
mass terms were not included. In that case, mixing occurs only amongst the n left-handed 
neutrinos that belong to the leptonic doublets, and there are one massive Dirac neutrino 
and (n — 1) massless neutrinos. Since a massive Dirac neutrino may be considered as 
a pair of Majorana neutrinos with identical mass, the model studied in Ref. [12] appears 
as a special case in the current work. The Majorana mass terms were indeed included 
in Ref. [11], but it was found that there are (2n — 1) mixing angles and n CP- violating 
phases in V c , totaling (3n — 1) free parameters, much more than the number found 
here. The difference arises from the fact that the constraints (22) from the texture zero 
in the original neutrino mass matrix have been completely exploited here to remove all 
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unphysical parameters while they were only partially applied in Ref.[12] to delete unitary 
transformations within the massless neutrinos. 

The above counting of physical parameters can be extended to the general case. With- 
out giving further details, we record the results as follows. For n > n > 1, V c contains 
n(no — 1) mixing angles and n(n — 1) CP phases, to be compared with n n + n(n — l)/2 
angles and n n + n(n — 1)/2 phases in Ref. [11]. Out of them, only n(2n — n — 1)/2 angles 
and n(2n — n— 1)/2 phases appear in V N . When n > n > 1, V° has n (n — 1) angles and 
n(n — 1) phases, much less than 2rion — no(no + l)/2 angles and 2n n — n — n (n — l)/2 
phases found in Ref. [11]. Out of those, only n (n — l)/2 angles and n (n — l)/2 phases 
enter in Vn- 

3 Phenomenological implications 

We study in this section some phenomenological implications of the leptonic CC interac- 
tions as revealed in eqs. (30,32,44). A possible way to measure the absolute neutrino 
mass is to study the electron spectrum in tritium (3 decay, which is sensitive to the effective 
neutrino mass 

"V = ^12 m i\ V v\ 2 = \V ln \^m 2 _c 2 m + m 2 + s 2 m = \V ln \y/m + m- = \V ln \\D\ (45) 

Thus the decay spectrum is only sensitive to Dirac-type mass in this model. 

The neutrinoless double (3 decay of nuclei is currently the only known practical means 
to unravel the Majorana nature of neutrinos. At the leading order of expansion in neu- 
trino mass over momentum transfer, the decay amplitude is proportional to the effective 
neutrino mass 

™W = lE m i(^5) 2 | = \Vm\ 2 \m-(c m ) 2 + m + {±is m f\ =0 (46) 
j 

Although the Majorana nature of neutrinos in principle allows the decay to occur, it is 
highly suppressed in the considered model. 

Finally, we want to study whether the model can accommodate the neutrino oscillation 
data excluding LSND [18]. This would require at least two mass squared differences and 
three mixing angles. From our above analysis, we know that with a single sterile neutrino 
there are only two free mixing angles for three generations. Although this is augmented by 
an effective mixing angle formed by the two neutrino masses, the chance to accommodate 
all data looks small. Our analysis below makes this claim stronger. As a first attempt, 
it is enough to consider the data that may be reasonably well described by vacuum 
neutrino oscillations, that is, those by KamLAND, K2K and CHOOZ. These experiments 
produce and detect neutrinos by CC interactions. The amplitude for the whole process 
that produces a neutrino and a charged lepton I a at the source and detects a charged 
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lepton ip at the detector is proportional to J2j Vaj*V/3j ex P[ — ^ m 'jL/(2E)], where L is the 
source-detector distance and E the energy of a relativistic neutrino. The probability is 

P(a ->/?)=*,*- 4 £ MQafin) sm 2 ^ + 2 ]T Im(Q Q/3;ij ) ^ ^ (47) 

i>j i>j 

where A - = m 2 - and Q a/3;lj = V^*yCyC* y c_ 

The above formula simplifies considerably in our model due to the special structure 
of V° and the existence of only three different masses. Eqs. (30, 32) give 

V £j V £* = E Wi, + K»V^- B + 4A n+ i) (48) 

7=1 

where the i factor in U drops out already. Taking into account the three different masses, 
we only need the following quantities 

n-1 



E <5a/3;n J = | V Qn | 2 (5 Q/3 - | Vs„ | 2 ) , 
3=1 

^2QaP;n+l,j = Sm\V an \ 2 (8 a p - \Vp n \ 2 ), 
3=1 

Qaf3;n+l,n = C m S m\Van,V(3n\ i 

where unitarity of V is used in the first two equalities. They are real even without using 
the reality of V in eq. (44). Thus the probability for the associated process with £ a , £p 
replaced by £ a , is the same. Denoting x± = m±L / (4E) , we obtain 

P(a->/3) = 5 a p-4\V an \ 2 (5 a p-\Vp n \ 2 )(cl l sm 2 x_ + s 2 m sm 2 x + ) . . 

-4:\V an V[3 n \ 2 c 2 m s 2 m sm 2 (x + - x_) 

Note that the probability does not normalize to unity [11]: 

£ P (a _ p) = 1 - 4 | V an \ 2 (c m s m ) 2 sin 2 (x + - x.) (51) 

This occurs because the mixing matrix appearing in the amplitude is not the one relating 
weak and mass eigenstates of neutrinos which is always unitary and would guarantee 
unity normalization, but the one appearing in CC interactions which is not unitary in the 
current model so that the sterile degree of freedom gets effectively lost in the sum. 

The KamLAND and CHOOZ are reactor v e disappearance experiments at L/E ~ 
(20 - 50) x 10 3 km/GeV and L/E ~ 333 km/GeV respectively, and are thus potentially 
sensitive to small and large mass squared differences. The K2K experiment observes 
accelerator disappearance at L/E ~ 200 km/GeV, close to the range at CHOOZ 
though in a different channel. That KamLAND and K2K observed deficits and spectral 
distortions implies that there would be two well-separated mass squared differences. Then, 
the null result by CHOOZ would be interpreted by some small mixing parameter. There 
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are two ways to arrange this, either (1) m 2 + ~ m 2 _ ^> rn^ — m 2 ^ or (2) m? + ^> m 2 _, implying 
correspondingly (1) c m ~ s m ~ 1/V2 or (2) c m ~ 1, s m ~ 0. In case (2), oscillations 
in the larger mass squared difference, i.e., in x + ~ x + — x_, are highly suppressed in all 
channels, thus incapable of accommodating the K2K data. In case (1), we have 

P(a -> a) ~ 1 - 4|\/ an | 2 (l - |\/ Qn | 2 ) sin 2 x + - |F Qn | 4 sin 2 (:r + - x_) (52) 

For K2K and CHOOZ, the last term can be ignored so that their results imply that 
4|V r 2n | 2 (l — \V 2 n\ 2 ) ~ 1 and 4 1 V" iri | 2 ( 1 — | Vl ri | 2 ) ~ 0. Since unitarity of V demands that 
\V\ n \ 2 + \V2n\ 2 < 1, the combined solution is, |V 2 „| 2 ~ 1/2 and \Vi n \ 2 ~ 0. But the latter 
is rejected by the KamLAND data. 

4 Conclusion 

We have investigated the leptonic mixing matrices in SM augmented by one sterile neu- 
trino. We found that the mixing matrix V c in CC interactions takes a factorized form, 
with one factor U describing the mixing in the subsystem of the two massive neutrinos 
and the other V the mixing of the massless neutrinos with the subsystem. We showed 
that this arises from the texture zero in the neutrino mass matrix that is protected by 
gauge symmetry. The matrix U is completely fixed by the two masses, while the fac- 
torization makes it possible to remove all phases in V, thus leaving us with (n — 1) 
free physical mixing angles in V° for n generations. The factorization also determines 
uniquely the mixing matrix V N in NC neutrino interactions in terms of the masses so 
that off-diagonal interactions occur only between the massive neutrinos. Both CC and 
NC interactions automatically conserve CP. We also considered some phenomenological 
results from the exposed structure in V c . The effective neutrino mass in tritium (3 decay 
is essentially sensitive to the Dirac mass in the model, while the leading contribution to 
neutrinoless double (3 decay vanishes. The difficulty with this simple model is that, even 
with any number of generations, it cannot accommodate the vacuum neutrino oscillation 
data coming from KamLAND, K2K and CHOOZ. A way out might be to introduce two 
sterile neutrinos, which would bring in more free parameters into the model. While the 
involvement of sterile neutrinos is not favored in solar and atmospheric data, it might still 
be phenomenologically viable by relaxing the tension between the LSND and other short 
baseline experiments, if either of the two can be unambiguously confirmed by MiniBooNE. 

References 

[1] For recent reviews on neutrino oscillations, see e.g.: M.C. Gonzalez-Garcia, Y. Nir, 
Rev. Mod. Phys. 75 (2003) 345; M. Maltoni, T. Schwetz, M.A. Tortola, J.W.F. Valle, 
New J. Phys.6 (2004) 122. 



12 



[2] For an introduction to neutrino physics, see e.g.: B. Kayser, F. Gibrat-Debu, F. 
Perrier, The Physics of Massive Neutrinos (World Scientific, 1989). 



[3] B. Pontecorvo, Sov. Phys. JETP 6 (1958) 429; Z. Maki, M. Nakagawa, S. Sakata, 
Prog. Theor. Phys. 28 (1962) 870; L. Wolfenstein, Phys. Rev. D 17 (1978) 2369; S.P 
Mikheyev, A.Yu. Smirnov, Sov. J. Nucl. Phys. 42 (1985) 913. 

[4] LSND Collaboration, C. Athanassopoulos et al, Phys. Rev. Lett. 75 (1995) 2650; C. 
Athanassopoulos et al, Phys. Rev. Lett. 77 (1996) 3082; C. Athanassopoulos et al, 
Phys. Rev. Lett. 81 (1998) 1774. 

[5] Bugey Collaboration, Achkar et al, Nucl. Phys. B 434 (1995) 503; CCFR Collabo- 
ration, Romosan et al, Phys. Rev. Lett. 78 (1997) 2912; KARMEN Collaboration, 
B. Armbruster et al, Phys. Rev. D 65 (2002) 112001. 

[6] BooNE Collaboration, E.D. Zimmerman et al, hep-ex/0211039. 

[7] M. Gell-Mann, P. Ramond, R. Slansky, in Supergravity (Ed., D. Freedman, P. van 
Nieuwenhuizen) (North- Holland, Amsterdam, 1979), p. 315; T. Yanagida, in Proceed- 
ings of the Workshop on Unified Theory and Baryon Number in the Universe (Ed. O. 
Sawada, A. Sugamoto) (KEK, Japan, 1979); R.N. Mohapatra, G. Senjanovic, Phys. 
Rev. Lett. 44 (1980) 912. 

[8] T.P Cheng, L.-F. Li, Phys. Rev. D 22 (1980) 2860. 

[9] J.T. Peltoniemi, D. Tommasini, J.W.F. Valle, Phys. Lett. B 298 (1993) 383; D.O. 
Caldwell, R.N. Mohapatra, Phys. Rev. D 48 (1993) 3259; Z.G. Berezhiani, R.N. 
Mohapatra, Phys. Rev. D 52 (1995) 6607; R. Foot, R.R. Volkas, Phys. Rev. D 52 

(1995) 6595; D. Suematsu, Phys. Lett. B 392 (1997) 413. 

[10] M. Maltoni, T. Schwetz, J.W.F. Valle, Phys. Rev. D 65 (2002) 093004, and its update 
cited in Ref. [1]. 

[11] J. Schechter, J.W.F. Valle, Phys. Rev. D 22 (1980) 2227. 

[12] J. Schechter, J.W.F. Valle, Phys. Rev. D 21 (1980) 309; see also, J.F. Donoghue, 
Phys. Rev. D 18 (1978) 1632. 

[13] A. Pilaftsis, Z. Phys. C 55 (1992) 275; B.A. Kniehl, A. Pilaftsis, Nucl. Phys. B 474 

(1996) 286. 

[14] KamLAND Collaboration, K. Eguchi et al, Phys. Rev. Lett. 90 (2003) 021802; T. 
Araki et al, hep-ex/0406035. 

[15] K2K Collaboration, M.H. Ahn et al, Phys. Rev. Lett. 90 (2003) 041801. 



13 



[16] CHOOZ Collaboration, M. Apollonio et al, Phys. Lett. B 420 (1998) 397; M. Apol- 
lonio et al, Phys. Lett. B 466 (1999) 415. 

[17] B. Kayser, A.S. Goldhaber, Phys. Rev. D 28 (1983) 2341; B. Kayser, Phys. Rev. D 30 
(1984) 1023; S.M. Bilenky, N.P Nedelcheva, S.T. Petcov, Nucl. Phys. B 247 (1984) 
61. 

[18] After submitting the paper to the arXiv, I was informed of the earlier work, F. del 
Aguila, J. Gluza, M. Zralek, Acta Phys. Polon. B 30 (1999) 3139, where the neutrino 
mass matrix was diagonalized explicitly for the case, n = 3, n = 1. The diagonalizing 
4x4 unitary matrix was used to compute the neutrino oscillation probability which 
thus normalizes to unity. However, as emphasized here and in Ref. [11], what really 
appears in the probability for a CC process is the rectangular 3x4 matrix in CC 
interactions instead of the above diagonalizing matrix even in the basis where the 
charged leptons are diagonalized. Nevertheless, their phenomenological conclusions 
for the case n — 3, n = 1 on the model's incapability of accommodating solar and 
atmospheric neutrino data are consistent with ours for any n. 



14 



